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Abstract. A novel method for calculating resonances in three-body Coulom- 
bic systems is presented. The Faddeev-Merkuriev integral equations are solved 
by applying the Coulomb-Sturmian separable expansion method. To show the 
power of the method we calculate resonances of the three-Q and the H~ sys- 
tems. 

1 Introduction 

For three-body systems the Faddeev equations are the fundamental equations. 
After one iteration they possess connected kernels, consequently they are effec- 
tively Fredholm integral equations of second kind. Thus the Fredholm alterna- 
tive applies: at certain energy either the homogeneous or the inhomogeneous 
equations have solutions. Three-body bound states correspond to the solutions 
of the homogeneous Faddeev equations at real energies, resonances, as usual in 
quantum mechanics, are related to their complex-energy solutions. 

The Faddeev equations were derived for short-range interactions and if we 
simply plug-in a Coulomb-like potential they become singular. The necessary 
modifications have been formulated in the Faddeev-Merkuriev theory [Q on a 
mathematically sound and elegant way via integral equations with connected 
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(compact) kernels and configuration space differential equations with asymp- 
totic boundary conditions. 

Recently, one of us has developed a novel method for treating the three- 
body Coulomb problem via solving the set of Faddeev-Noble and Lippmann- 
Schwinger integral equations in Coulomb-Sturmian-space representation. The 
method was elaborated first for bound-state problems Q with repulsive 
Coulomb plus nuclear potential, then it was extended for calculating p — d 
scattering at energies below the breakup threshold ||. Also atomic bound- 
state problems with attractive Coulomb interactions were considered [Q. In 
these calculations an excellent agreement with the results of other well estab- 
lished methods were found and the efficiency and the accuracy of the method 
were demonstrated. 

More recently we have extended the method for calculating resonances in 
three-body systems with short-range plus repulsive Coulomb interactions by 
solving the Faddeev-Noble integral equations ||. Here we solve the Faddeev- 
Merkuriev integral equations. This way we can handle all kind of Coulomb-like 
potentials, not only repulsive but also attractive ones. For illustrating the power 
of this method we show our previous three-a results and, as novel feature, we 
calculate resonances in the H~ (pee) system. 

2 Faddeev-Merkuriev integral equations 

The Hamiltonian of a three-body Coulombic system reads 

H = H° + W + v a +Vp+v 7 , (1) 

where H° is the three-body kinetic energy operator, W stands for the possible 
three-body potential and v a denotes the Coulomb-like interaction in the sub- 
system a. We use throughout the usual configuration-space Jacobi coordinates 
x a and y a . Thus v a only depends on x a (v a = v a (x a )), while W depends on 
both x a and y a coordinates (W = W(x a ,y a )). 

The physical role of a Coulomb-like potential is twofold. Its long-distance 
part modifies the asymptotic motion, while its short-range part strongly corre- 
lates the two-body subsystems. Merkuriev proposed to split the potentials into 
short-range and long-range parts in the three-body configuration space via a 
cut-off function £, 

v^(x a) y a ) = v a (x a )(; a (x a ,y a ), (2) 

and 

v < £>(x ai y a ) = v a (x a )[l - C, a {x a ,y a )]. (3) 

The function £ a is defined such that it separates the asymptotic two-body sector 
fl a from the rest of the three-body configuration space. On the region of fl a 
the splitting function £ Q asymptotically tends to 1 and on the complementary 
asymptotic region of the configuration space it tends to 0. Rigorously, fl a is 
defined as a part of the three-body configuration space where the condition 

\x a \ < x (l + \y a \/y ) 1/u ,with x > 0, y > 0, v > 2, (4) 
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is satisfied. So, in Q a the short-range part v a s coincides with the original 
Coulomb-like potential v a and in the complementary region vanishes, whereas 
the opposite holds true for v a . Note that for repulsive Coulomb interactions 
one can also adopt Noble's approach H, where the splitting is performed in the 
two-body configuration space. This approach can be considered as the yo — > oo 
limit of Merkuriev's splitting. Then v a coincides with the whole Coulomb 

(s) 

interaction and v a with the short-range nuclear potential. 

In the Faddeev procedure we split the wave function into three components 

W) = Ha) + Hp) + Ha), (5) 

where the components are defined by 

H a ) = G^(z)v^\V). (6) 

Here G"' is the resolvent of the long-ranged Hamiltonian 

ffW =H° + W + vV) +v { < ) +vV\ (7) 

G^(z) = (z — H^ l ')~ l , and z is the complex energy parameter. The wave- 
function components satisfy the homogeneous Faddeev-Merkuriev integral 
equations 

H a ) = G^(z)v^J2\^ ( 8 ) 

for a = 1, 2, 3, where Gq' is the resolvent of the channel long-ranged Hamilto- 
nian 

H {i) =fr (0 +u W ! ( 9 ) 

Ga (z) = (z — Ha'*)' 1 . Merkuriev has proved that Eqs. (||) possess compact 
kernels for positive E energies, and this property remains valid also for complex 
energies z = E — iT/2, r > 0. 



3 Solution method 



We solve these integral equations by using the Coulomb-Sturmian separable 
expansion approach. The Coulomb-Sturmian (CS) functions are defined by 



(r|n) 



(n + 2l + l)\ 



1/2 



(2br) l+1 exp(-br)L% +1 (2br), 



(10) 



with n and I being the radial and orbital angular momentum quantum numbers, 
respectively, and b is the size parameter of the basis. The CS functions {\n)} 
form a biorthonormal discrete basis in the radial two-body Hilbert space; the 
biorthogonal partner defined by (r\n) = r _1 (r|n). Since the three-body Hilbert 
space is a direct product of two-body Hilbert spaces an appropriate basis can 
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be defined as the angular momentum coupled direct product of the two-body 
bases (the possible other quantum numbers are implicitly assumed) 

\nv) a = \n) a ® \v) a , {n,v = 0,1,2,...). (11) 

With this basis the completeness relation takes the form 

N 

1 = lim Y, »<H = lim l£. (12) 

n,v=Q 

Note that in the three-body Hilbert space, three equivalent bases belonging to 
fragmentation a, (3 and 7 are possible. 

In Ref. ||] a novel approximation scheme has been proposed to the Faddeev- 
type integral equations 

|^) = GW(^iM s) E 1 tK>» (13) 

i.e. the short-range potential v a in the three-body Hilbert space is taken to 
have a separable form, viz. 

JV 

where = a (n^|wi s ' \n'v')p. The validity of the approximation relies on the 

square integrable property of the term u„ |i/; 7 ), 7^a. Thus this approxima- 
tion is justified also for complex energies as long as this property remains valid. 
In Eq. ( |l4| ) the ket and bra states are defined for different fragmentation, de- 
pending on the environment of the potential operators in the equations. Now, 
with this approximation, the solution of the homogeneous Faddeev-Merkuriev 
equations turns into solution of matrix equations for the component vector 

^jj_ a = a (m>\lp a ) 

± a =Q^{z)v^Y±^ ( 15 ) 

where Gffi = a (jw\Ga \n'v') a . A unique solution exists if and only if 

det{[G ( ' ) (z)]" 1 -i; (s) } = 0. (16) 

The Green's operator Gh} is a solution of the auxiliary three-body prob- 
lem with the Hamiltonian H& ■ To determine it uniquely one should start again 
from Faddeev-type integral equations, which does not seem to lead any further, 
or from the triad of Lippmann-Schwinger equations Q . The triad of Lippmann- 
Schwinger equations, although they do not possess compact kernels, also define 
the solution in an unique way. They are, in fact, related to the adjoint repre- 
sentation of the Faddeev operator The Hamiltonian H a l \ however, has a 
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peculiar property that it supports bound state only in the subsystem a, and 
thus it has only one kind of asymptotic channel, the a channel. For such a sys- 
tem one single Lippmann-Schwinger equation is sufficient for an unique solution 

I- 

The appropriate equation takes the form 

Gg> ^G a + G a U a G^, (17) 

where G a is the resolvent of the channel-distorted long-range Hamiltonian, 

H a = H +v a + ug), (18) 

and U a = W + Vp + — uf . The auxiliary potential depends on the 
coordinate y a and has the asymptotic form u„' ~ e a (ep + e 7 )/y Q as y a — > oo. 
In fact, is an effective Coulomb interaction between the center of mass of 
the subsystem a (with charge ep + e 7 ) and the third particle (with charge e a ). 
Its role is to compensate the Coulomb tail of the potentials lA + u 7 in fl a . 

It is important to realize that in this approach to get the solution only the 
matrix elements G^ are needed, i.e. only the representation of the Green's 
operator on a compact subset of the Hilbert space are required. So, although 



Eq. (17) does not possess a compact kernel on the whole three-body Hilbert 
space its matrix form is effectively a compact equation on the subspace spanned 
by finite number of CS functions. Thus we can perform an approximation, 



similar to Eq. (|lj), on the potential U a in Eq. (17), with bases of the same 
fragmentation a applied on both sides of the operator. Now the integral equa- 
tion reduces to an analogous set of linear algebraic equation with the operators 
replaced by their matrix representations. The solution is given by 

[G^izT^U^z^-ir. (19) 

The most crucial point in this procedure is the calculation of the matrix 
elements G a — a (nv\G a \n'v') ai since the potential matrix elements v *l and 
U_ a can always be calculated numerically by making use of the transformation 
of Jacobi coordinates. The Green's operator G a is a resolvent of the sum of 
two commuting Hamiltonians, H a = h Xa + h Va , where h Xa — h x + v a and 

h Va = hy a + Ua , which act in different two-body Hilbert spaces. Thus, us- 
ing the convolution theorem the three-body Green's operator G a equates to a 
convolution integral of two-body Green's operators, i.e. 

G a {z) = -L I dz'g Xa {z-z') g ya {z'), (20) 
J c 

where g Xa {z) = (z — fr^) -1 and g Va (z) = (z — hy^ 1 . The contour C should 
be taken counterclockwise around the continuous spectrum of h Va so that g Xa 
is analytic in the domain encircled by C . 

To examine the structure of the integrand let us shift the spectrum of g Xa 
by taking z = E + is with positive e. By doing so, the two spectra become 
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well separated and the spectrum of g Va can be encircled. Next the contour C 
is deformed analytically in such a way that the upper part descends to the 
unphysical Riemann sheet of g Va , while the lower part of C can be detoured 
away from the cut [see Fig. |l|. The contour still encircles the branch cut sin- 
gularity of g Va , but in the e — > limit it now avoids the singularities of g Xa . 
Moreover, by continuing to negative values of e, in order that we can calculate 
resonances, the branch cut and pole singularities of g Xa move onto the second 
Riemann sheet of g Va and, at the same time, the branch cut of g Va moves onto 
the second Riemann sheet of g Xa . Thus, the mathematical conditions for the 
contour integral representation of G a (z) in Eq. (^) can be fulfilled also for 
complex energies with negative imaginary part. In this respect there is only a 
gradual difference between the bound- and resonant-state calculations. Now, 
the matrix elements G Q can be cast in the form 

G a (z) = ±-l dz'g x Jz - z>) g ya (z% (21) 
J c 

where the corresponding CS matrix elements of the two-body Green's operators 
in the integrand are known analytically for all complex energies |^] , and thus 
the convolution integral can be performed also in practice. 



4 Numerical illustration 



4.1 Resonances in a model three-alpha system 

To show the power of this method we examine the convergence of the results 
for three-body resonant-state energies. For this purpose we take first the same 
model that has been presented by us before in Ref . || . This is an Ali-Bodmer- 
type model for the charged three-a system interacting via s-wave short-range 
interaction. To improve its properties we add a phenomenological three-body 
potential. Adopting Noble's splitting we have 

„W(r) = V x expl-r 2 //^ 2 } + V 2 exp{-r 2 //3 2 2 } (22) 

with V x = 125 MeV, V 2 = -30.18 MeV, fa = 1.53 fm, fa = 2.85 fm, and 

v «(r) = 4e 2 /r. (23) 

We use units such that h 2 /m = 41.47 MeV, e 2 = 1.44 MeV fm. The mass of the 
a-particle is chosen as M — 3.973m, where m denotes the mass of the nucleon. 
The three body potential is taken to have Gaussian form 

W(p) =Fexp{V//3 2 }, (24) 

3 

where p 2 = J2 v h V = -31.935 MeV and j3 = 3.315 fm. Here r, stands for the 

i=X 

position vector of i-th particle in the center of mass frame of the three-a system. 
Since we consider here three identical particles we can reduce the necessary 
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Faddeev components to one. We select states with total angular momentum 
L = 0. In Table 1 we show the convergence of the energy of the ground state 
and of the first resonant state with respect to N, the number of CS functions 
employed in the expansion. The selected resonance is the experimentally well- 
known sharp state which has a great relevance in nuclear synthesis. 

4.2 The H~ system 

In this system all the two-body interactions are of pure Coulomb type, two of 
them are attractive, therefore we have to use the genuine Merkuriev approach. 
Since the two electrons are identical particles we can reduce the number of 
Faddeev components to two. Furthermore, we can define the cutoff function 

(s) 

such that Ve e ' = 0. In this case the corresponding Faddeev component vanishes 
identically. So, finally we have to deal with one Faddeev component only. In 
the pe subsystem we take the Merkuriev's cut with the functional form 

C(x, y) = 2/(1 + exp((a:/a:o)7(l + v/Vo))), (25) 

and with the actual parameters are v — 2.1, xq = 5 and yo = 10. Fig. 2 and 3 
show the short and long range parts of v pe , respectively. In Table 2 we present 
the convergence of the energy of a L — resonant state with respect to N and 
the number of angular momentum channels used in the bipolar expansion. 

5 Conclusions 

In this article we have presented a new method for calculating resonances in 
three-body Coulombic systems. The homogeneous Faddeev-Merkuriev integral 
equations were solved for complex energies. For this, being an integral equation 
approach, no boundary conditions are needed. We solve the integral equations 
by using the Coulomb-Sturmian separable expansion technique. The method 
works equally for three-body systems with repulsive and attractive Coulomb 
interactions. 
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Figure 1. Analytic structure of g Xa [z — z') g Va (V) as a function of z with z = E + ie, 
E > 0, e > 0. The contour C encircles the continuous spectrum of h Va . A part of it, 
which goes on the unphysical Riemann-sheet of g Va , is drawn by broken line. 



Table 1. Convergence of the ground-state and of the first resonant-state energy (in 
MeV) of a three-a system interacting via the potential of (^2|) with increasing basis for 
the separable expansion. N denotes the maximum number of basis states employed 
for n and v in Eq. (0). 
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Table 2. Convergence of a resonant state of the H~ system with respect to the 
number of CS states. The angular momentum channels are taken into account up to 
I = 2, I = 3 and I — 4, respectively. The (complex) energies are measured in atomic 
units. 
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